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Abstract: We consider the standard Abelian sandpile process on the Bcthc 
lattice. We show the existence of the thermodynamic limit for the finite volume 
stationary measures and the existence of a unique infinite volume Markov process 
exhibiting features of self-organized criticalityf 1 1. 



1 Introduction 

Global Markov processes for spatially extended systems have been around for 
about 30 years now and interacting particle systems have become a branch of 
probability theory with an increasing number of connections with the natural 
and human sciences. While standard techniques and general results have been 
collected in a number of books such as [Liggett (1985)[ |Chen (1992)1 , |Toom (19901 



and are capable to treat the infinite volume construction for stochastic systems 
with locally interacting components, some of the most elementary questions for 
long range and nonlocal dynamics have remained wide open. We have in mind 
the class of stochastic interacting systems that during the last decade have in- 
vaded the soft condensed matter literature and are sometimes placed under the 
common denominator of self-organizing systems. 

Since the appearance of the paper ]BTW (1988) ], the concept of self-organized 



criticality (SOC) has suscitcd much interest, and is applied in a great variety of 
domains (see e.g. fTurcotte (1999)| for an overview). From the mathematical 
point of view, the situation is however quite unsatisfactory. The models ex- 
hibiting SOC are in general very boundary condition dependent (especially the 
BTW model in dimension 2), which suggests that the definition of an infinite 
volume dynamics poses a serious problem. Even the existence of a (unique) 
thermodynamic limit of the finite volume stationary measure is not clear. From 
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the point of view of interacting particle systems no standard theorems are at 
our disposal. The infinite volume processes we are looking for will be non-Feller 
and cannot be constructed by monotonicity arguments as in the case of the 
one-dimensional BTW model (see [MRSV (2000) |) or the long-range exclusion 
process (see [ Liggett (19801 ). On the other hand in order to make mathemati- 
cally exact statements about SOC, it is necessary to have some kind of infinite 
volume limit, both for statics and for dynamics. 

In this paper we continue our study of the BTW-model for the case of the 
Bethe lattice, this is the abelian sandpile m odel on an in finite tree. For this 
system, many exact results were obtained in | DM (1990b) |. In contrast to the 
one-dimensional case this system has a non-trivial stationary measure. We show 
here that the finite volume stationary measures converge to a unique measure [i 
which is not Dirac and exhibits all the properties of a SOC-state. We then turn 
to the construction of a stationary Markov process starting from this measure 
^. The main difficulty to overcome is the strong non-locality: Adding a grain at 
some lattice site x can influence the configuration far from x. In fact the cluster 
of sites influenced by adding at some fixed site has to be thought of as a critical 
percolation cluster which is almost surely finite but not of integrable size. The 
process we construct is intuitively described as follows: At each site x of the 
Bethe lattice we have an exponential clock which rings at rate <p(x). At the 
ringing of the clock we add a grain at x. Depending on the addition rate <p(x), 
we show existence of a stationary Markov process which corresponds to this 
description. We also extend this stationary dynamics to initial configurations 
which are typical for a measure fi' that is stochastically below /i. 

The paper is organized as follows. In section 2 we introduce standard results 
on finite volume abelian sandpile models, and summarize some specific results 
of pM (1990b)] for the Bethe lattice which we need for the infinite volume 
construction. In section 3 we present the results on the thermodynamic limit 
of the finite volume stationary measures and on the existence of infinite vol- 
ume Markovian dynamics. Section 4 is devoted to proofs and contains some 
additional remarks. 



2 Finite Volume Abelian Sandpiles 

In this section we collect some results on abelian sandpiles on finite graphs which 
we will need later on. Most of these results are contained in the review paper 
[pliar (1999)1, or in pP (1998)|. 



2.1 Toppling Matrix 

Let V denote a finite set of sites and = (A x y ) x , ye v a matrix indexed by 
the elements of V satisfying the conditions : 

1. For all x, y £ V, x + y, A^ y = A^ < 0, 

2. For all x 6 V, A^ x > 1, 

3. For all x e V, E ye v ^ly > 0, 
4- L,^^>0. 
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Such a matrix is called a toppling matrix. The fourth condition ensures that 
there are sites (so called dissipative sites) for which the inequality in the third 
condition is strict. This is fundamental for having a well defined toppling rule 
later on. In the rest of the paper we will choose A v to be the lattice Laplacian 
with open boundary conditions on a finite simply-connected set V C S, where 
S is a regular graph, like the d-dimensional lattice Z d , or the infinite rootless 
treeTd of degree d+ 1. More explicitly: 

A^ x = 2d if V C 1 d } 

= d+liiVcT d , 
A x = — 1 if x and y are nearest neighbors. (2.1) 

The dissipative sites then correspond to the boundary sites of V. 



2.2 Configurations 

A height configuration 77 is a mapping from V to IN = {1,2, ...} assigning to 
each site a natural number r/(x) > 1 ("the number of sand grains" at site x). A 
configuration 77 £ ¥S V is called stable if, for all x 6 V, rj(x) < Ajf x . Otherwise 
rj is unstable. We denote by fly the set of all stable height configurations. For 
n £ TN V and V C V, rj\y denotes the restriction of 77 to V. 

2.3 Toppling Rule 

The toppling rule corresponding to the toppling matrix A v is the mapping 



defined by 



T A v(n,x)(y) = n(y) - A^ y if rj(x) > A^ x , 

— n(y) otherwise. (2.2) 

In words, site x topples if and only if its height is strictly larger than A x x , by 
transferring ~A x y grains to site y ^ x and losing itself A x x grains. Toppling 
rules commute on unstable configurations. This means for x, z G V and 77 such 

that v (x) >Al x , n{z) >A% tZ , 

T A v (T A v(r],x),z) =T A v (T A v(r),z),x) (2.3) 

We write [T A v (•, z)T A v (•, x)](n) = [T A v (•, x)T A v (•, z)\{rf). 

Choose some enumeration {x\, • • • , x n } of the set V. The toppling transfor- 
mation is the mapping 

defined by 



T A v : m v Q v 



T AV ( V ) = lim Y[T A v(; Xl ) (77). (2.4) 



In ]IP (1998)| it is recalled that 
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1. The limit in (2.4) exists, i.e. there are no infinite cycles, due to the presence 
of dissipative sites. 

2. The stable configuration T&v {rj) is independent of the chosen enumeration 
of V . This is the abelian property and follows from (2.3). 



2.4 Addition Operators 

For 77 6 TN V and x <G V, let rj x denote the configuration obtained from r\ by 
adding one grain to site x, i.e. rj x {y) — n(y) + 5 XjV . The addition operator 
defined by 

a x ,v '■ — > £V; V >-> a x yr] = T A v (rj x ) (2.5) 
represents the effect of adding a grain to the stable configuration rj and letting 



the system topple until a new stable configuration is obtained. By (2.3), the 



composition of addition operators is commutative: For all n € f2y, x,y G V, 

a x ,v{a v yn) = a y ,v{a x>v r)). 



2.5 Finite Volume Dynamics 

Let p denote a non-degenerate probability measure on V, i.e. numbers p x , 
< p x < 1 with JjjjgyPa = 1. We define a discrete time Markov chain 
\ji n : n > 0} on f2y by picking a point x € V" according to p at each discrete 
time step and applying the addition operator a x .v to the configuration. This 
Markov chain has the transition operator 

Pvf(v) = Y,P*K a *> vr l)- ( 2 - 6 ) 

x£V 

We can equally define a continuous time Markov process {r/ t : t > 0} with 
infinitesimal generator 

Kfiv) = E vWt/K,^) - /(»?)], (2.7) 

x£V 

generating a pure jump process on with addition rate tp(x) > at site x. 



2.6 Recurrent Configurations, Stationary Measure 

We see here that the Markov chain {n n ,n > 0} has only one recurrent class and 
its stationary measure is the uniform measure on that class. 

Let us call IZy the set of recurrent configurations for {r] n , n > 0}, i.e. those 
for which P v (n n = rj infinitely often) = 1, where P v denotes the distribution of 
{r) n ,n > 0} starting from r/o = r\ G fiy. In the following proposition we list 
some properties of Hy. For the sake of completeness we include a proof which 
we could not find worked out completely in the literature. 

Proposition 2.1 

1. TZy contains only one recurrent class. 

2. The composition of the addition operators a x y restricted to TZy defines 
an abelian group G. 
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3. \G\ = \K V \. 

4- For any x G V , there exists n x such that for any rj G TZy , a ™V^ = r l- 
5. \K V \ = det A v . 

Proof: 1. We write r\ ( if in the Markov chain £ can be reached from rj with 
positive probability. Since sand is added with positive probability on all sites 
(p x > 0), the maximal configuration rj max defined by 

Vmax( x ) = A-x.x 

can be reached from any other configuration. Hence, if 77 € TZ y then ?y <^-> rj max , 
therefore r\ max G TZy and ?7 maa ; <^-> 77 (see e.g. [ Chung (1960) p. 19). 
2. Fix 77 £ 7£y; then there exist n y > 1 such that 

n 

and 



.'/< ",.r 11 "„.\ 

satisfies {a x yg x )(ri) = (g x a>x,v)(v) = V- Tne set 

={(£^: (tfe.vfcXC) = C} 

is closed under the action of a x ,Vi contains T), hence also n max : it is a recurrent 
class. By part 1, 1Z X = TZy, a x yg x is the neutral element e, and g x = a~\r if 
we restrict a x y to TZy. 

3. Fix ( e Ry and put * c : G -> TZ v ;g h-> g(Q. As before * C (G) is a 
recurrent class, hence ^(G) = TZy. If for g, h G G, ^(.g) = ^(/i)) then 
gh^ 1 ^) — C, and by commutativity gh~ 1 (g'() — g'C, for any g 1 G G. Therefore 
gh^ 1 ^) — £ for all £ G TZy, thus g = h. This proves that ^ is a bijection from 
G to fty. 

4. Since G is a finite group, for any x £ V there exists n x > 1 such that 
",V = '■ 

5. Adding A^ x particles at a site x G V makes the site topple, and — A x ^ y 
particles are transferred to y. This gives 



_ TT 

x,V - 11 a y.V 

y^x 



On TZy the a x y can be inverted and we obtain the closure relation 

11 a S v = '• 

y£V 

which completely determines the one-dimensional representations of the group 
of addition operators, and in particular the cardinality of the latter, as obtained 



in jDhar (1990a 



Remark. TZy does not depend on th e p x , an d do es not change by going from 
discrete to continuous time, i.e. from ( |2.6| ) to (2.7). 

The main consequence of the group property of G is the fact that the unique 
stationary measure is uniform on TZy. 
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Proposition 2.2 

1. The measure 

is invariant under the action of a x y, x S V (6 V is the Dirac measure on 
configuration rj). 

2. On L 2 (jiy) the adjoint of a x ,v * s 

a l,v = a xy- ( 2 - 9 ) 
Proof: Since a x y ■ Tiv Hv can be inverted, we have 



hence (2.9). By choosing g = 1, part 1 follows. ■ 

Remark. This shows that fiy is invariant under the Markov processes gener- 
ated by (O) and (O). 



2.7 Burning Algorithm 

The burning algorithm determines whether a stable configuration n £ Sly is 
recurrent or not. It is described as follows: Pick n e ily and erase all sites 
x G V satisfying the inequality 

rj(x)> £ (- A D- 

This means "erase the set Ei of all sites x £ V with a height strictly larger than 
the number of neighbors of that site in V" . Iterate this procedure for the new 
volume V\Ei, and the new matrix A v \ El defined by 

A^ = Al y ifx,yeV\E 1 
— otherwise, 

and so on. If at the end some non-empty subset Vf is left, rj satisfies, for all 
x 6 Vf, 

The restriction r)\y f is called a forbidden subconGguration. If Vf is empty, the 
configuration is called allowed, and the set Ay of allowed configurations satisfies 

Proposition 2.3 



Ay = TZy. 
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The main ingredient to prove this result (see [ IP (1998)] , |Spccr (1993) |) is 



the fact that toppling or adding cannot create a forbidden subconfiguration. 
The set Ay is thus closed under the dynamics and contains the maximal con- 
figuration r) max - 

Remark that the burning algorithm implies that for V' D V, and 77 G Oy, 
77 £ TZy, then any ( G Sly such that C|v = V satisfies £ ^ 72. y Indeed, 
the property of having a forbidden subconfiguration in Vf only depends on 
the heights at sites x G Vf. Therefore r\ G 72.y implies 77 |y G 72y. This 
"consistency" property will enable us to define allowed configurations on infinite 
sets. 

2.8 Expected Toppling Numbers 

For x, y G V and 77 G f2y, let y, 77) denote the number of topplings at site 

7/ G V by adding a grain at x G V, i.e. the number of times we have to apply 
the operator T&v(-,y) to relax rf . Define 



G v (x,y) = J fiv(dTj) n v (x, y, 77). (2.10) 
Writing down balance between inflow and outflow at site y, one obtains (cf. 



[|Dhar (1990a)| ) 

A x,zGv(z,y) = S x , y , 

z£V 

which yields 

G v (x,y) = (A v )-* y . 

In the limit V f S (where S is 7L d or the infinite tree), Gy converges to the 
Green's function of the simple random walk on 5*. 

2.9 Some specific results for the tree 

When V n is a binary tree of n generations, many explicit results have been 



obtained in | DM (1990b) |. We summarize here the results we need for the 



construction in infinite volume. 

1. When adding a grain on a particular site G V n of height 3, the set of 
toppled sites is the connected cluster (73(0,77) of sites including having 
height 3. This cluster is distributed as a random animal (i.e. its distribu- 
tion only depends on its cardinality, not on its form). Moreover 

lim n Vn QC 3 (0,V)\ = fc) - Ck- 3 / 2 (2.11) 

n\ 00 

as k goes to infinity. The notation ~ means that if we multiply the left 
hand side of ( |2.11 ) by fc 3 / 2 , then the limit k — » 00 is some strictly positive 



constant C. 



2. When adding a grain on site x, the expected number of topplings at site 
y satisfies 

i™ 1 / Mv„(*7) n Vn (x,y,rj) = G{x,y), (2.12) 

"Too J 
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where G(x, y) is the Green's function of the simple random walk on the 
infinite tree, i.e. 

G(0,x) =C2~W, (2.13) 
and \x\ is the "generation number" of x in the tree. 

3. The correlations in the finite volume measures \xy h can be estimated in 
terms of the eigenvalues of a product of transfer matrices. This formalism 
is explained in detail in [DM (1990b) |, section 5: Let /, g be two local func- 
tions whose dependence sets (see below a precise definition) are separated 
by n generations. To estimate the truncated correlation function 

/A4 {f;g) = J fgdn Vk - J fd^ Vk J gd^ Vk , (2.14) 

consider the product of matrices 

n / i i k,n -i i k,n \ 

^=n( 1+ r ("« 

where 7^'™ € [0,1]. Let AJ^ fc (resp. A^} fc ) denote the smallest (resp. 
largest) eigenvalue of M%. Then 

i*vAf;9)<c(f,g)^s- (2-16) 

If / and g have dependence sets "deep within" v\, then 7. is very close 
to one, and the correlations are governed by the maximal and minimal 
2 2 



eigenvalues of M n = 



1 3 



3 Main results 

3.1 Notation, definitions 

From now on, S denotes the infinite rootless binary tree, l^cSa finite subset 
of S; fly is the set of stable configurations in V, i.e. fly = {t] : V — > {1, 2, 3}}, 
and the set of all infinite volume stable configurations is £1 = {1,2,3} S . The 
set Q is endowed with the product topology, making it into a compact metric 
space. For 77 6 il, n\v is its restriction to V, and for 77, ( 6 f2, 77y(y c denotes the 
configuration whose restriction to V (resp. V c ) coincides with r/\v (resp. C|i /C )- 
As in the previous section, IZy C fV is the set of all allowed (or recurrent) 
configurations in V , and we define 

K = {nen-.W C S finite, r/\ v G ^y}. (3.17) 

A function / : — ► 1R is Jocai if there is a finite V C S such that 77 |y = (\v 
implies f(rj) = /(C)- The minimal (in the sense of set ordering) such V is called 
dependence set of / and is denoted by Df. A local function can be seen as a 
function on Qy for all V D Df and every function on Qy can be seen as a local 
function on f2. The set £ of all local functions is uniformly dense in the set C(fi) 
of all continuous functions on f2. 

All along the paper, we use the following notion of limit by inclusion for a 
function / on the finite subsets of the tree with values in a metric space (K, d) : 
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Definition 3.1 Let S = {V C 5, V finite}, and f : S -> (isT,d). Then 

lim f(V) = k 

if for all e > 0, i/iere exists Vq 6 S suc/i i/ia£ /or a/Z D Vb, d(f(V), k) < e. 



Definition 3.2 ^4 collection of probability measures vy on fly is a Cauchy net 
if for any local f and for any e > there exists Vq D Df such that for any 
V,V'DV 

f(v)vv(dn) - J f{rj)vv'{dri)\ < e. 

A Cauchy net converges to a probability measure v in the following sense: 
The mapping 

* : C -» R; / *(/) - lim / /(fry 

defines a continuous linear functional on C (hence on C(f2)) which satisfies 
> for / > and = 1. Thus by Riesz representation theorem 

there exists a unique probability measure on such that ^f(f) = J fdv. We 
denote vy — > 2/, and call this ^ the infinite volume limit of iy- 

We will also often consider an enumeration of the tree S, {xq, x\, . . . , x ni . . .}, 
and put 

T n = {x , . . . ,x n }. (3.18) 
3.2 Thermodynamic limit of stationary measures 



Theorem 3.1 The set 1Z defined in (3.17) is an uncountable perfect set, i.e. 

1. 1Z is compact, 

2. The interior oflZ is empty, 

3. For all rj G 1Z there exists a sequence Tj n ^ r\, n n 6 7Z, converging to n. 

For T) G we denote by (73(0,77) the nearest neighbor connected cluster of 
sites containing the origin and having height 3. 



Theorem 3.2 The finite volume stationary measures fly defined in (2.1 ) form 
a Cauchy net. Their infinite volume limit fi satisfies 

1. fi(TZ) = I, 

2. (i is translation invariant and exponentially mixing, 

3. /x(r ? :|C 3 (0,r7)|<oo) = l } 

4. J\C 3 (0,v)\Kdv) =00. 

Remark: Point 3. above remains true for the set Ci(0, 77), the nearest neighbor 
connected cluster of sites containing the origin and having height 1, and probably 
also for 62(0,77) but this we have not been able to prove. 
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3.3 Infinite volume dynamics 

The finite volume addition operators a x .v (cf. ( |2.5| )) can be extended to ft via 
a x y : ft — * ft; n i-> a^yr? = (a KiV -^|v)y??y c ■ (3.19) 

Proposition 3.1 

1. There exists a subset ft' ofTZ with n(fl') = 1 on which the limit 

lim a x vV = a x r\ (3.20) 
vts ' 

exists, and a x r\ G ft! . 

2. The measure p of theorem \3.2\ is invariant under the action of a x . 

3. For every r\ G ft', a x (a y r]) = a y (a x rf), for all x,y G S. 



Part 2 implies that the infinite volume addition operators a x (cf. ( 3.20| )) 
define norm 1 operators on L p (/i), for 1 < p < 00 via 

(axf)(v) = f{axV)- 

We now construct a Markov process on /it-typical infinite volume configura- 
tions which can be described intuitively as follows. Let ip : S — > (0, 00); this 
function will be the addition rate function. To each site x G S we associate 
a Poisson proces N^ x (for different sites these Poisson processes are mutually 
independent) with rate f{x). At the event times of N^ x we "add a grain" at x, 
i.e. we apply the addition operator a x to the configuration. Then Ly introduced 
in generates a pure jump Markov process on ft. Indee d, t his operator is 



well-defined and bounded on any L p (n) space by Proposition 3J, which implies 



Proposition 3.2 L v is the L p {fi) generator of the stationary Markov process 
defined by 



exp(tl%)f = J \T[a" r f) dP 



where P denotes the joint distribution of the independent Poisson processes 
N^,andfeLP(p). 

The following condition on the addition rate ip is crucial in our construction. 
Remember \x\ is the generation number of x: 



Summability Condition: ^ <p(x)2 I 35 ' < 00 (3.21) 



xes 



This condition ensures that the number of topplings at any site x G S 
remains finite during the addition process. 



Theorem 3.3 If ip satisfies condition (3.21), then we have 
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1. The semigroups Sy(t) — exp(iiy) converge strongly in L (/J>) to a semi- 
group Sip(t). 

2. S v (t) is the L 1 ^) semigroup of a stationary Markov process {r]t : t > 0} 
on f2. 



3. For any f € C, 



t|0 t 



where the limit is taken in L 1 ^). 



Remarks. 

1. In Proposition 4.1, we show that S v (t) is a strongly continuous function of 
<P- 

2. In Proposition 4.2, we show that condition (3.21) is in some sense optimal. 



Theorem 3.4 The process {rj t : t > 0} of Theorem 3.S admits a cadlag version 
(right- continuous with left limits). 

The intuitive description of the process {r/ t : t > 0} is actually correct 
under condition (3.21), i.e. the process has a representation in terms of Poisson 
processes: 

Theorem 3.5 If ip satisfies condition ( 3.2\ ), for fixW almost every (r],uj) the 
limit 



lim I 

i/t.<? J- J- 



0.1 



V = rjt 



x£V 



exists. The process {r]t '■ t > 0} is a version of the process of Theorem 3.S, i.e. 
its i 1 (/i) semigroup coincides with S v (t). 



Finally, we can slightly generalize Theorem 3.5 in order to define the dy- 
namics starting from a measure stochastically below /i. For 77, £ £ il, 77 < £ 
if for all x G S, rj(x) < C( x )- A function / : O — > JR is monotone if ?/ < ( 
implies f{rf) < /(C)- Two probability measures /1 and v satisfy fx < v if for all 
monotone functions, J fd[i < J fdv. 



Theorem 3.6 Let // < ji. If tp satisfies condition ( 3.21 ). for // x P almost 
every {r],uj) the limit 



lim TT 

i/ts J- J- 



V1S 



x£V 



exists. The process {i] t : t > 0} is Markovian with rjo distributed according to 



Remark. The last Theorem implies that 77 = 1 can be taken as initial configu- 
ration. 
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4 Proofs 

This section is devoted to the proofs of the results described above. Some of 
them will be put in a slightly more general framework so that they can be ap- 
plied to other cases (where S is not a binary tree or where we have other addition 
operators a x ) as soon as the existence of a thermodynamic limit of the finite 
volume stationary measures is guaranteed. The essential cause of difficulty is 
the non- locality of the addition operators. The essential simplification is the 
abelian property which enables us to think of the a x as complex numbers of 
modulus one. 



4.1 Thermodynamic limit of stationary measures 



Proof of Theorem |3.1| : 

1,2. If rj G TZ and ( > r], then ( £ TZ (by the burning algorithm £|y > rj\v 
implies that £|y G TZy). Since rj = 2 is in TZ (again by the burning algorithm), we 
conclude that TZ is uncountable. To see that TZ has empty interior, notice that if 
rj G TZ, there does not exist x, y G S nearest neighbors such that rj{x) = i](y) = 1 
(that way, ri\[ Xty y would be a forbidden subconfiguration). Finally TZ is closed 
as intersection of closed sets. 

3. Let rj max be the maximal configuration, r/ maa; (a;) = 3 for all x G S. If 
r/\v G TZy, then rjy (r] max )v <= G TZ. Therefore any rj G TZ containing aninfi- 



nite number of sites x for which r/(x) ^ 3 has property 3 of Theorem 3.1. If 
?/ G TZ contains only a finite number of sites having height 1 or 2, then we 
choose a sequence E = {x n : n G IN} C {x G S : r/(x) — 3 and rj{y) = 
3 for any neighbor of x} such that two elements of £ are never nearest neigh- 
bors, and \x n \ is strictly increasing in n. We then define r] n (x) = rj(x) for 
x G S \ {xk G £ : < k < n} and r] n (xk) = 2 for xu G E, k > n+ 1. These r\ n 
belong to TZ by the burning algorithm, and r\ n — > rj. ■ 



Proof of Theorem 3.2 



We use T n introduced in ( 3.18| ), but with the Xi such that n < m implies that 
the generation numbers satisfy | I . Then we have 



log 2 n. 



(4.22) 



To prove that the probability measures /xy form a Cauchy net, it is sufficient to 
show that for any local function / : — > IR we have 



fdfiT n ~ fdfJ-T n+1 \ < oo. 



(4.23) 



We do it for f(rf) = f](xo) (a general local function can be treated in the same 
way), by giving an upper bound of the difference J fd[iT n — J fd^T n+1 by a 



truncated correlation function (cf. (2T4J)). Then we estimate the latter by the 
transfer matrix method (cf. Section 2.9, part 3). We abbreviate in what follows 
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Lemma 4.1 

iMn+ifafco)] -Mn[^(^o)]| < C/x n+ i[7?(a;o);I(7?(x n+ i) = 3)] 

Proof: By the burning algorithm, every r\ £ 7?. n can be extended to an element 
of 7£„+i by putting rj(x n +i) — 3. Moreover 



thus 

which yields 



{vIt,, ■ V € 7?.„ + i, r?(a;„ + i) = 3} = Tl n , 



fi n+ i[r](x. 



n+lj 



3] 



n+l 



(4.24) 



E Jl^o) 



E 
E 



ri(x )I(ri(x n+1 = 3) 



+H 



n+l\ 



-r](x )I(rj(x 



ra+lj 



3)- 



Therefore 



iMn+iMzo)) - Mn(??(a:o))| < 



Mu+lh^n+l) = 3] ' 

fJ-n+i[y{xo); I(rj(x n+1 ) = 3)] 
H n+1 [r]{x n + 1) = 3] 



The Lemma follows now from ( 4.24 ), and the fact that |7?. ra | grows like e cn for 
some c > log 2. ■ 



Recalling Section 2.9, part 3, we have 



^n[n{xo);I{r](x n ) = 3)] < C 



X„\,\x n \ 



l*„l,k«l 



A 



(4.25) 



M 



Lemma 4.2 



+ oo .|a;„|,|x n | 
Am 



E 



|x n | ,\x n | 



< +oo 



Proof: We abbreviate Am' 1 = Am" 

^|x„|,|a„|^ j^ emem ber < ^j^™''' 1 ™ < 1 is close to one for i <C n and n large. In 
terms of the trace and the determinant of M(n) we have 



^•■'*»',M(n) = M£ n|)7 i 



\(") 
A (") 



i (Tr(M(n)) + \/[Tr(M(n))] 2 - 4det(M(n)) 
i (Tr(M(n)) - \/[Tr(M(n))] 2 - 4det(M(n)) 



Therefore, 



/A^\ / [Tr(M(n))] 2 
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To prove the Lemma we show that (cf. (fO|)) 

/ det(M(n)) \ 4 |8w| 
vjTr(M(n))]V " V ' 



Use 



and 



det(M(n)) = nii+^O 2 



! = 1 



Tr(M(n)) > Tr J| 



l + 7z 

2 + 7l ; 



n( i+ ^)+n( 2+ ^)' 



to estimate (for 1 < i < \x n \, 2(2 + 7,) > 3(1 + 7;)) 



(4.26) 



det(M(n)) 
[Tr(M(n))] 2 



< 



i +2 n 



2 + 7, 



l + 7< 



n 



2 + 7. 
l + 7i 



< (l + 2.(3/2) l "" l + (3/2) 2|a: " 1 )" 1 

< (4/9) 1 *" 1 . 

For a general local function /, we have to replace \x n \ by \x n \ — No, where 
TVo is the number of generations involved in the dependence set of /. Since / is 
local, iVo is finite, hence the convergence in ( 4.23| ) is unaffected. ■ 



4.2 Infinite volume toppling operators 

Definition 4.1 Given the finite volume addition operators a x .v (defined in 



(3.H)) acting on f2, we call a configuration r\ £ Q normal if for every x G S 



there exists a minimal finite set V x (j]) C S such that for all V D V D V x (tj) 

a x ,V'V = a x,W- 

In other words, for a normal r\, outside V x (j)), no sites are affected by adding 
a grain at x. In our case, when a particle is added at some site x £ S, the 
cluster of toppled sites coincides with the cluster Cs(x, rj) of sites having height 
3 including x, thus 

V x (r l )=C 3 (x,r l )Ud e C 3 (x,r 1 ), (4.27) 

where d e denotes the exterior boundary. Notice that for a normal configuration 
77, by definition, 

Qx(v) = ^o, x ,v(v) = a x ,v x (r,){v) ( 4 - 2 8) 
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Proof of Proposition 3.1 



1. We show that there is a full measure set Q' of normal configurations. From 
( p.lip and Theorem ^72 , 



I(\C 3 (x,r,)\ =n)rf^C^ 2 . 



Therefore concentrates on the set f2' of configurations for which all the clusters 
C 3 (x,rj) are finite, hence for which r\ is normal. Moreover this set is closed 
under the action of the addition operators a y , since (cf. (4.27)) 



C 3 (x,a v r}) cV x (r))UV y (r)), 



(4.29) 



2. Choose e > 0, pick a local function /, fix V n f S and no such that n > no 
implies 



l-i{r] £ f2 : V x (rj) t V n } < 



4II/IU + 1 

This no exists since /i concentrates on normal configurations. We estimate 
f(a x r])dfj, - /(n)d/i| < 



(4.30) 



f(a x ,v n v) d ^ - J .f(v)d/J.\ 
211/Hoo/ifo e n : T4(r?) £ V n } 



< lim 

m 

e 



< 77 + 2||/||oohm^y m (a X! v„ ))) ^ a^y^n)) 

= | + 2||/||oo (l - linw m (W C K)) 

= | + 2||/|| 0O (l-/i(K(n)cK))<e. 

In the last step we used that the indicator I(V X (rj) C V n ) is a local function. 

3. Let 7] 6 O', x,y £ S be two different sites and V Z) V x (n) U V x (a x r]) U 
V x (a y rj). Since a x y and a y- y commute, we have 

a x (a y rj) = a x (a y yn) = a x y(a y yrj) 
= a v ,v(ax,vrj) = a y.v( a xV) = a y (a x rj). 



4.3 Infinite volume semigroup 



We now turn to the proofs of Theorems 3.2 and 3.4. 
Definition 4.2 We define the cluster of rj £ at x E S as 
C(x,rj) = {yeS: a y v(x) ^ rj(x)}, 

and put 

Gn(x,y)= / I(y £C(x,r/))dfj,(r]). 



(4.31) 
(4.32) 
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Finally for ip : S — > [0, oo), write 

WfWv = ^2 l fi( x ) / V>{dr{)\f{a x rj) - f{r)% 
xes J 

B v = {/ : f2 — ► IR : / bounded, \\f\\ v < oo}. 

Lemma 4.3 J/ 

ip(x)G fl (y,x) < oo for all y E S, (4.33) 
then all local functions are in B v . 

Proof: Let / be a local function with dependence set Df. Then f(a x rj) ^ f{rj) 
if for y G D f , a x n(y) ^ i](y), i.e. xGC(y,rj): 

ll/llp = ^2<P( X ) / \ a xf - fW 

= J vi x )Wf - fW 

x£U yeDf C(y,r)) 

xes J 
< 211/Hoo ^ £G„(i/ j!B MaO<oo. 



The next lemma provides a link between and the Green's function for 



simple random walk on S, i.e., between conditions (4.33) and (3.21). 
Lemma 4.4 

G^(x, y) < Y G(y, z) = 6 XiV + 3G(x, y), 

z^x 

where z ~ x means that z and x are neighbors. 

Proof: We have to estimate the probability that a x n(y) ^ r)(y). If by adding 
a grain at x we influence y, this can only be achieved by the toppling of one of 
the nearest neighbor sites of y. Since fi concentrates on normal configurations, 

m ( a xv(y) v(y)) = ^^(^(y) ¥= v(y),v x (j]) uv y (rj) c v) 

= Jim lim (a x ,w(y) ¥= V(v), v x(v) U V y {rj) C V) 
= lim lim (*w (a x ,wv{y) ¥= v(v), V x (r]) U V y (n) C V) 

v \ j W \ o 

< ^ s ^w(a x ,wv(y)¥=v(y)) 

< Urn n w (3 z G W, z - y,n w (z,y,r)) > 1) 

- IPoE / d ^w(v)nw(z,y,v) 

^G(z,y), (4.34) 



Z~j, 
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where we used (|TTot) , (pTTl) , Cgg) and flOgl), 



The following Lemma finishes the proof of Theorem 3.3 and shows that £> v is 
a natural core for the domain of the generator of the infinite volume semigroup. 

Lemma 4.5 

1. For f £ B v the net 



S*{t)f = exp(iL£)/ = exp It £ - I) j / 



(4.35) 



converges in (as V ] S) to a function S v (t)f S / i— > S v (t)f 

defines a semigroup on B v which is a contraction in both L 1 ^) and B v 
norms. 



2. Under condition ({.3c), the semigroup S v (t) corresponds to a unique Markov 
process on Q. 

Proof: We denote by ||/|| the L 1 ^) -norm of /, and we abbreviate Sv{t) = 
S£(t),S(t) = Sv(t),L v = L*. 

1. First note that Sv(t) is well-defined on by Proposition 3.2. By the 

abelian property (Proposition [Tl|, part 3) we can write for V C V C S: 

\\S v (t)f - Sv>W\\ = \\(Sy\v(t) - I)Sv(t)f\\ 



By Proposition 3.2, Sv(t) is the semigroup of a stationary Markov process and 
hence a contraction on £ (/u). Therefore 

\\S V (t)(Sy,\ V (t) - I)f\\ < \\(S V ,\ V (t) - I)f\\ 

= II / L V '\vSy\v(s)fds\\ 
Jo 

< f \\L v ,\ v f\\d8 

Jo 

< t J2 <P( X ) /lK -I)f\dn-+0 as V,V T S, 



(4.36) 



xev\v 



where the last step follows from / 6 B v . Hence Sv{t)f — > S(t)f in L 1 ^). We 
show that S(t)f £ B v : 



xes ^ 
< ^ p(s) f S(t)|ox/-/|d/x 

= J^2<p(x)\a x f-f\d^=\\f\\ v> . 



Thus is also a contraction for the || • || y -norm. We finish with the semigroup 
property: 

S(t)S(s)f = Ijm S v (t)[S(s)f] 

= $Tsms Sv{t)Sw{t)L 
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and 



S(t + s)f = KmS v (t)S v (s)f. 



Then, since Sv(t) is a contraction in L^/i), 

\\S v (t)S w (t)f - S v (t)S v (s)f\\ < \\S w (s)f - S v (s)f\\, (4.37) 



By ( }4.36| ) , the ri ght h and side of ( |4.37| ) goes to zero as V, W | S 1 . 
2. If condition ( 1.33] ) is met, then B v contains all local functions by Lemma 
1.3. Therefore, by contractivity the semigroup S(t) on B v uniquely extends 
to a semigroup of contractions on L x (/i). Since by Proposition 3.2, <SV(i) is a 
Markov semigroup, so is S(t), i.e. S(t)l = 1, S(t)f > if / > 0. Hence, by 
Kolmogorov's theorem there is a unique Markov process with semigroup S(t). ■ 



Remark. When ip = l, condition ( 1.35 ) is equivalent to 

X J K dr i)I(^ e C(y,rj)) = J \C(y,r))\(i(dr)) < +oo, 



i.e., the clusters must be integrable under /i. For models which exhibit "self- 
organized criticality" , C(y, rj) is usually a "finite but critical percolation cluster" , 
implying that / \C(y, rj)\dn = oo (cf. Theorem 2^2 part 4, because C(y,n) D 
deC^^y,!])). Therefore this formalism breaks down for addition rate (p = 1. 



The following Lemma proves Theorem 3.4 



Lemma 4.6 Under condition ( ^.3!\ ), the process {r) t : t > 0} o/ Theorem 
is almost surely right- continuous, i.e. 



P„ 



limd(ry t ,ry ) > e 



0. 



(4.38) 



where P^ is its path-space measure, and the distance d is defined below (in 



Proof: Pick a function : S — > (0, 1) such that 



xES 



and 



The distance 



<P(x)G»(x,v)*(v) 



< oo 



d(r ) ,Q = '£\r 1 (x)-C(x)\nx) 



(4.39) 
(4.40) 
(4.41) 



x£S 



generates the product topology. Denote by E M the expectation w.r.t. P M . For 

fy(v) = v(y)> t 

fv(Vt) - fy(vo) = [ L*f y ( Vs )ds + Mf, 



(4.42) 



where is a centered martingale with quadratic variation 



[(M? ) 2 ] = IE, 



(L*tf( Vs )-2f y (r ls )L*f y ( Vs ))ds 



Abclian Sandpile Model 



19 



Using stationarity of rj s and 



dii\L*g\ <2\\g\\ 00 J2 £ p(s)G„(y, 

x6S ye£> s 



for a local bounded function on f2, we obtain from ( 4.42 ) 



E„ [(Mf) 2 ] <C^^)G^(j/ )a :). 



Now we can estimate 



P, 



3s<i:^|?/ s (y)-r ?0 (y)|*(y)>e 

yeS 

/ d S y;|L*'/ 1 ,(» ?a )|*(!/)>e/2 
70 yes 

< (12t/e) ]T rtx)G^y,x)*(y) + (2/e) 2 TE 



< P 



sup 

0<s<t 



x, yes 



yes 

2 



> e/2 



yes 



< (12i/e) £ V (x)G fl (y,x)^(y) + (2/e) 2 TE, £(M» ) 2 *( 2/ ) 

x,yeS yes 

a^yeS 

Here we used Markov's and Doob's inequalities in the second step and the 



Cauchy-Schwarz inequality combined with (4.3E) in the third step. The result 
(pi follows. ■ 



4.4 Poisson representation 



In this section we prove Theorems p.5| and 3.6. Intuitively it is clear from the 



abelian property that the process of which we showed existence in the previ- 

ous subsection can be represented as rLeS Qx * Vi where N^ x are independent 
Poisson processes of intensity <p(x). 

We take T n as in ( 3.18 ). We say that the product Y\ xeS a^"q exists if for 
every y £ S there exists N y such that for all m, n > N y 



n <^ 



(y) 



n <^ 

,i6T m 



(y) 



o. 



This is equivalent to the convergence of the sequence II^eT a x XT ) m the product 
topology. 

Lemma 4.7 Under condition fy.3!\ ), the product 

i [ ax* V = fjt 

xes 
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exists for /i-almost every realization of Ny x and almost every rj. The process 
{fit ■ t > 0} is a version of the Markov process of Lemma 4-' ■ 

Proof: Choose a realization of N*' x such that 



(4.43) 



xes 



for every y. This happens with probability one by condition ( 4.33| ). Define for 
r) e CI' 



Vrjt)= [1 a. 



(4.44) 



xET„ 



Under fi, rjT n (t) is stationary in n and t. We have 

n[\(vT n (t))(y)-(vT n+ M)(y)\ >i] 

(tf)- far. (*))(») 



< 



a x * +1 VT n (t) 



li{dri) 



N v " 



[y) - v(y) 



< 



/ Yj I [ a L+^] (v) 

J 3=1 



7 — 1 



(y) /i(djj) 



n+lj 



In the second and last steps we used the invariance of \x under a x . By the 
Borel Cantelli Lemma, ( [4.43 ) implies that for almost every realization of N^ x 



H [3n : Vn > n {r)T n (t))(y) = (VT no (t))(y)] 



1. 



This proves /i-a.s. convergence of the product. To see that fjt is a version of the 
Ma rkov process with semigroup S v (t), combine Proposition 3^ with Theorem 
3.3, part 1 to get, for any local function /, 



dft 



dPf(f} t )-S v (t)f 



= 0. 



In the preceding argument we used a particular enumeration of the countable 
set S. But changing it gives again a process with semigroup S v (t). Therefore 
the limiting process will not depend (up to sets of measure zero) on the chosen 
enumeration of S. ■ 



Proof of Theorem 3.6: 



For r) 6 O' and y G S we have the relation (remember ( 4.44[ ) ) 

w(t)(y) = v (y) + i(y eV)Y, N T ~ AnV(v), 



(4.45) 



where n t v (x) 1 an integer valued random variable, is the number of topplings at 
site x in the time interval [0, t], when sand is added in V. For T n defined in (3.18) 
we will first prove that increases /*xP almost surely to an L 1 (jtixP) random 
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variable n*, interpreted as the number of topplings in [0, t] when we add grains 
according to N^ x . By the abelian property the sequence n^ fc (0) is increasing in 
k. The following estimate is similar to (4.34) 



(p X F) ( K fc (0) - n^ +1 (0) | > e) = ( M X F) (0) > e) 

n* (0)^(d»j)xP(tL;) 



< 



1 



< ~ limy < +1 (0)/(^ +1 Wuy M c V)n{dn) x p(dw) 

< - lim / n\ v (x k+1) 0, rf)p w (drj) x P(dw) 

< -^(z fe+1 )G(0, Xfc+i). (4.46) 
e 

In the fourth line, n^(a;fe+i , 0, rf) denotes the number of topplings up to time t 
at site S W by adding grains at site irfc+i £ W. By the Borel Cantelli Lemma, 
condition ( [3.21 ) implies the a.s. convergence of n^ fe (0), and analogously of every 
n T k ( x )- Pick (?7, w) such that n l Tk (r), u>) converges, i.e. such that sup fc n l Tk (77, u>) — 
n*(77, w) is hnite (indeed, rij, (77, uj) is an integer). If 77' < 77, then rij. (r]',ui) < 
n l Th (77, bS) because we can obtain r\ from 77' by adding sand at sites ieS for which 
T)'(x) < r)(x) thereby increasing the number of topplings. We thus conclude that 
the convergence of (77, w) implies the convergence of n^ k (77', oj) for all 77' < 77. 

Now let p! < p in the FKG sense. There is a coupling F12 of p! x F and 
/j x F such that 

F12 (((7/1, ojx), (772,^2)) : u>i = ll>2,Vi < m) = ! 5 

i.e. we use the same Poisson events and couple p' and p according to the optimal 
coupling (see fStrassen (1965) |). Then 

(// x F) (n^(i),w) -Wfaw)) = P i2 (n^fa.wi) -m'fe^)) 

= Pl2 (nT k (Vl, ^l) ^*('7l ? '*'i),7lr k (j72,W2) -» 7l t (?72, W 2 ), CJi = W 2 ,?7l < ^2) 

> F12 (n^fe) <*> 2 ) -> n t (?72,a;2)) 

= (^xP)(n^(»7,w)->n*(t7 > w)) = l. 



This shows the // x F-almost sure convergence of n^ fc , hence by ( 4.45 ) the prod- 
uct Yl x&s dx v ^ 77' converges p! x F almost surely. ■ 

As a further result we show that the semigroup S v (t) is continuous as a 
function of the addition rate <p. We define 

i x = {if : S -> [0, 00) : \\ip\\ = ^ <p(x)G(0, x) < 00}. 

It is a complete metric space (as a closed subset of a Banach space) with the 
property: If <p n G ^1, Pn T <P (pointwise), and ip € £1 then y>„ — > ip in £1. 



Proposition 4.1 The semigroup S v (t) of Theorem 3.c is a strongly continuous 
function of ip, i.e. if tp n — > ip in l\, then for any local function f, S Vn (t)f — > 

s v {t)f. 
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Proof: Let 



firm 



be the number of topplings in [0, t] from sand 



addition at rate ip. In the proof of Theorem 3.5 we have shown that this random 
variable is /i x P almost surely well defined and, after taking limits in Q4.45| ), 
satisfies 



= m + N'-Ar 



where JV£ 



limy TS J^xev N v' X - Note that if - 



(4.47) 

the Poisson processes 



and can be coupled in such a way that for all x G S, Nj^ < Nj*, and 
hence, by the abelian property, riL^x) < n^ 2 (x). Consider a coupling of the 



1 N* AVn ,N*v Vn **dN* 
Xi(t) > X 2 {t), X 2 (t) < X 3 (t), X 3 (t) > X 4 (t), are satisfied with probability one. 



four Poisson processes iV*, N vAif>n , ]V! V „ and iV* under which the inequalities 



Let P d enote the law of the marginal (Xi,X±) 

to (pi), 



We have, by a reasoning similar 



dji(E|n* Vn (0)-n$,(Q)|) < / dfi (e (n^„(0) - nJ, nA „(0)) 

+ E« Vv „(0)-< AVji (0)) 

+ E(n*, VpB (0)-nJ,(0))) 

< ^|^„(x)-^(x)|G(0,x), (4.48) 



xGS 



which tends to zero for tp n — ► tp in l\. Take now a local function /, and denote 
D f =D f Ud e Df, 



\S Vn (t)(f) - S v (t)(f)\ < W(n% n {x)^nl{x) 



for some x £ D 



< nn%{x)-n% n {x) 



x£D, 



Combining this with ( 4.48j ) concludes the proof. 



One might ask whether we can go beyond condition fl3.2l|) , which essentially 
guarantees that the expected number of topplings stays finite in the addition 
process. In the following proposition we show that it is impossible to keep 
integrable toppling numbers and "rate 1" a dditi on. The relation (4.49|) should 



be regarded as the infinitesimal version of (4.47), where a(x) replaces the rate 
<p(x). We then show that <p has to depend on x. 

Proposition 4.2 Let a : S — > {0,1} be a stationary and ergodic random field 
distributed according to v. Denote by J a{Q)v{drj) — p its density. Suppose 
there exists a measurable transformation T : {0, 1} S x Q — > f2 which satisfies the 
conditions 



1. The measure /i of Theorem S.i is invariant under T(a, •) for any a. 
2. 

T(a, rf)(x) = r/(x) + a(x) — An(a, rj, x), 
with n(a,r], .) 6 L 1 ^) for v almost every a. 



(4.49) 



Then, p = 0. 
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Proof: Taking expectation over p in (4.49|) gives 



a(x), 



(4.50) 



where ^{a,x) = J n(a,T),x)p,(drf). By stationarity of p and i/, ty(a,x) is a 
stationary random field. Let (x t : t > 0) denote continuous time simple random 



Divide this last line by t and let t j +oo. As ^ is ergodic (making the last term 
equal to p) and as the process ^(a, xt) is stationary, we conclude that p = 0. ■ 
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